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Abstract 

The initial states which minimize the predictability loss for a damped harmonic oscillator 
are identified as quasi-free states with a symmetry dictated by the environment's diffusion 
coefficients. For an isotropic diffusion in phase space, coherent states (or mixtures of coherent 
states) are selected as the most stable ones. 



PACS numbers: 03.65.Bz, 05.30.-d, 05.40.+j 

The emergence of classical reality from the underlying quantum description of the world is 
one of the most fascinating unsolved problems of present-day physics. Decoherence was proposed 
as a mechanism for the selection of classical (preferred) states: due to the interaction with an 
environment (external degrees of freedom) classical states are singled out as the most stable ones. 

Zurek, Habib and Paz |] invented a criterion, called "predictability sieve", for distinguishing 
the preferred set of states from the rest of the Hilbert space: classical states are characterized by 
the least increase in entropy. They addressed this problem in the context of the Caldeira-Leggett 
model (a particle moving in a potential and linearly coupled with a bath of harmonic oscillators, 
|], 0, 0, @ ) • They found out that coherent states are selected, via predictability sieve, as the most 
"classical" ones. 

In contradistinction to this approach, we will study the problem of classicality in the framework 
of Lindblad's theory, in which the structure of the master equation is derived from very general 
assumptions concerning the mathematical description of the evolution of an open quantum system 
||. This strategy leads to general results (since Lindblad equations work for a large class of 
physically interesting systems |J), and also avoids any problems related to the non-preservation 
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of the positivity of the reduced density matrix, which was perceived by some authors [0] as an 
inconvenient of the Caldeira-Leggett type of evolution. We will identify the states with classical 
behavior as the states which minimize the rate of entropy increase. We will apply this criterion in 
an identical fashion for both pure and mixed initial states, using a formalism that doesn't require 
one to discriminate between them. 

Lindblad's result shows that a general form for the generator of a completely positive dynamical 
semigroup in the Schrodinger picture is || 

Hp) = -{\h,p\ + |r£([^ + ] + Wj,pv?]), (l) 

3 

where for the operators Vj, we will take || Vj — ajp + bjq, where j = 1,2 while a,j and bj 
are complex numbers. For the Hamiltonian H of the system we will consider a second-order 
polynomial operator in coordinate and position written as a sum of a harmonic-oscillator part and 
a mixed term 

muj 2 9 p 2 a , 

We denote the diffusion coefficients by 

h 2 h 2 Ti 2 

D qg=2^2 N 2 ' D pp = 2 ^ ' fei ' 2 ' Dpq = Dqp = ~2 Re ^ a * jbj ' ^ 
j=i j=i j=i 

and the friction constant is 

A = -Im ][>*&,, (4) 
i=i 

It is easy to show |J that the following inequality must be satisfied 

D pp D qq -D 2 pq >^-. (5) 

With the above notation, the master equation which governs the evolution of the system takes 
the form 

^ = ~[H,p(t)]-j([q,p(t)p]-\p,p(t) q ]) 

-^[p, \p,p(t))) - [q,p(t))) + ™fZ\p, [q,p(t))). (6) 

For the correlations between two operators C and C we will use the definition 

*c,c>(t) = Tr [pW \ j - <rc(t)a c >(t), (7) 

where acif) — Tr(p(t)C) and ac'if) — Tr(p(t)C). will denote the dispersion matrix 
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With the notations above, Heisenberg's inequality reads detH(t) > For £(£) the time- 

evolution is known |j) 

^ = yS(t) + S(t)F T + 2I?, (9) 

where 

-(A — /i) a; 



Y 



-uj — (A + /x) 



r is the transposed matrix of 1" and T> is a 2 x 2 matrix with elements = muD qq , T> pq = 
T^qp = D qp and £> pp = 

For the case of the damped harmonic oscillator, it is known that a certain class of states 
(quasi- free states ||) has the property of invariance under the action of the quantum semigroup. 
This can be seen by writing the Fokker-Planck equation which corresponds to (H) in the form || 

df w (x 1 ,x 2 ,t) * d 1 2 d 2 

St = Y ij — {x j f w {x Xl x 2 ,t)) + — ^ V ij ——f w {x 1 ,x 2 ,t), (10) 

where fw is the Wigner function of a quasi-free state 

fw(xi,x 2 , t) = [(27r) 2 detS(t)]-^ exp |— ^X T (t)E _1 X(i)| , (11) 

and X(t) = ( Xl a<1 /l \ It is now easy to verify that the Gaussian Wigner function (|TT]), with 

the time-dependence of the mean values cr p (£), cr q {t) and of dispersions o~ qq (t), o~ pp (t) , <J pq (t) given 
by equations (3.26 - 3.27) from 0, is a solution of the equation (pi]). Thus, quasi-free states are 
preserved during the evolution of the system. 

Starting with such a state, we are interested to calculate the rate of linear entropy variation. 
The initial states with the lowest rate of linear entropy increase will be identified as the most 
stable (i.e. the most classical-like) states. 

The linear entropy is a convenient measure of the purity of a quantum state and is defined by 

S (p) = l-Tr(p 2 ). (12) 

For a quasi-free state we have 

s( P ) = l-j. (13) 



where A(t) is the "area" in phase space, A(t) = jr^detTi(t). The condition for a quasi-free state 

to be a pure one is A{t) = 1, that is detY*{t) = ^ (the equality case in the Heisenberg relation). 
The time-derivative of s(t) can be calculated by making use of the following relation 

d " n&ffl "» = Tr f^E-'( ( )) , (14) 

where 

IpVf) = - | ~maT ~ CX pq( t ) \ 

w detE{t) \-a pq (t) muja m {t) ) ' v ; 



Then we obtain from (Q) and (]I3J) that 

dA(t) 2 



dt h 

But the rate of the linear entropy increase is given by 

ds(t) 1 dA(t) 
dt ~ A{t) 2 dt 



detY>{t){TrY + Tr{VT,- l {t))). (16) 



(17) 



so the behavior of the rate of linear entropy increase is given entirely by the time derivative (|16|) 
of the area A(t) for all states (including pure states). In the following we shall find the squeezing 
parameter of the initial state for which the rate of increase of A{t) is minimized. 

Let us notice, for the beginning, that the positive 2x2 matrix S(t) can be diagonalized, at 
any particular instant t, in the form (see also [T2|, |T3| |) 

hA(t) T ftt 2 (t) 



where K(t) is a real positive number (the squeezing parameter) A(t) is the area occupied by the 
system in phase space and 0(t) is an orthogonal symplectic matrix for which we will employ the 
usual form 

(t)=( cos6 { t ) ~ sin9it A 



sin0(t) cos0(£) J 
A similar formula holds for the D-matrix 



2 u V d 

with 



P = ^OSff ° )0 D , (19) 



-» f (2Q) 



sin (p cos (p 

Here, A is a parameter controlling the intensity of diffusion, d characterizes the degree of anisotropy 
and ip is the rotation angle. 
Now, QTBD and QT7D imply 



ds 
di 



A(Q) 



' 2A + Tr(S- 1 (0)P) . (21) 



This result is even more general when the system is in a pure state at t=0 (so that A(0) = 1), 
in the sense that it does not depend on the kind of initial states we are starting with (quasi-free 
states or not). Indeed, in |J it was shown that 

= 2Tr{pL{p)) = ~ J2(Tr(pV jP V;) - Tr(p 2 V;VA) > 0. (22) 
For a pure state p 2 = p and pOp = Tr(pO)p for any selfadjoint operator O. Then 

= I Ed^WOI 2 - TripVfV,)) > 0. (23) 
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We have Tr(pVj) = a^a p + bj<r q . Hence we get 



-2A + -^(D qq a pp {0) + D qq a pp {0) - 2D pq a pq {0)) 



-2A + Tr [Tj (0)V 



which has the same form as the rate of linear entropy increase calculated before (see ([H]) ) for pure 
(A(0) = 1) initial quasi- free states. 

We are interested in finding the states which produce the least increase of the phase-space area 
at the initial moment t = 0. By minimizing the expression 



— 2A + Tr (E _1 (0)X> 



-2A 



A 

W) 



{cos^(0) - ip)[K 2 (0)d- 2 + X- 2 (0)d 



+ sin 2 (#(0) - v)[K 2 (0)d 2 + *r 2 (0)<T 2 ]}, 



with respect to K(0) and 0(0), one gets 



nun 



ds 
~dt 



corresponding to 



t = 



N*(0) 
0*(O) 



A- A{0)X 

! A{oy ' 



: d, 
tp. 



(24) 



(25) 



So, in the general case of an anisotropic diffusion, the minimum variation of the area in phase 
space is obtained when the squeezing parameter of the state equals d, the degree of anisotropy 
of the diffusion, and the characteristic rotation angle of the correlation and diffusion matrices are 
equal. For an isotropic diffusion, V pp = V qq , V pq = 0, i.e. d=l (and the rotation angle vanishes 
trivially from all the relations, since now the system has rotation symmetry in phase space); we 
obtain K*(0) = 1. This case corresponds to many models of dissipation, especially from quantum 
optics (see the same result was obtained by Zurek, Habib and Paz |TJ in the context of 

the Caldeira-Leggett environment model. In other words, a phase-space isotropic environment 
favors a symmetric state, while an anisotropic environment selects states with the same degree of 
anisotropy. 

When the initial state is pure, 



nun 



ds 
~d~t 



t = 



2 (A - A) > 0, 



(26) 



where the last inequality comes from (|5|) and expresses the fact that the linear entropy of pure 
states in an environment always increases, because the pure initial states become more and more 
mixed. 

Our result shows that the values of K*(0) and 6*(0) are independent of the overall magnitude 
A of the diffusion. They are also independent of A(0). Thus, the same degree of squeezing is 
singled out, irrespective of the purity of the initial state, thus confirming previous insights 0, |J 
regarding the structure of the mixed preferred states: they can be seen as the thermalization of 
the selected pure states. 
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We conclude by emphasizing the main results of this paper. In general, the pure or mixed state 
which produces the minimum rate of increase in the area occupied by the system in phase-space 
is a quasi- free state which has the same symmetry as that induced on the evolution in phase-space 
by the diffusion coefficients. For isotropic phase-space diffusion, the selected pure states are the 
coherent states. 
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